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Abstract 

Electron accelerators with higher and higher longitudinal field gradients are desir- 
able, as they allow for the production of high energy beams by means of compact 
and cheap setups. The new laser-plasma acceleration technique appears to con- 
stitute the more promising breakthrough in this direction, delivering unprecedent 
field gradients up to TV/m. In this article we give a quantitative description of the 
impact of longitudinal wake fields on the electron beam. Our paper is based on 
the solution of Maxwell's equations for the longitudinal field. Our conclusions are 
valid when the acceleration distance is much smaller than the the overtaking length, 
that is the length that electrons travel as a light signal from the tail of the bunch 
overtakes the head of the bunch. This condition is well verified for laser-plasma 
devices. We calculate a closed expression for the impedance and the wake function 
that may be evaluated numerically. It is shown that the rate of energy loss in the 
bunch due to radiative interaction is equal to the energy emitted through coherent 
radiation in the far-zone. Furthermore, an expression is found for the asymptotic 
limit of a large distance of the electron beam from the accelerator compared with the 
overtaking length. Such expression allows us to calculate analytical solutions for a 
Gaussian transverse and longitudinal bunch shape. Finally, we study the feasibility 
of Table-Top Free-Electron Lasers in the Vacuum Ultra- Violet (TT-VUV FEL) and 
X-ray range (TT-XFEL), respectively based on 100 MeV and 1 GeV laser-plasma 
accelerator drivers. Numerical estimations presented in this paper indicate that 
the effects of the time-dependent energy change induced by the longitudinal wake 
pose a serious threat to the operation of these devices. 
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1 Introduction 



The quest for the ultimate electron accelerator will never end. A better beam 
quality, a higher energy, a smaller size and lower costs will always be goals 
to be pursued. In particular, size and costs are strictly related. This justifies 
research activities towards higher and higher longitudinal field gradients. 
The availability of higher longitudinal field gradients is obviously related, 
in its turn, to the possibility of production of beams with higher energies. 

At the present day, particle accelerator technology is mainly based on Radio- 
Frequency (RF) devices. However, in the past few years a novel laser-plasma 
acceleration method has been developed and experimentally demonstrated 
IH B @]/ that promises to outclass all existing accelerator technologies. The 
laser-plasma acceleration technique requires an ultrashort (a few femtosec- 
onds long) high intensity (Tera to Peta-Watt) laser pulse focused into a 
supersonic Helium gas jet (or filled capillary) with density around 10 19 /cra 3 . 
This produces a region free of electrons propagating behind the laser pulse. 
Many electrons (up to 10 10 ) are captured in this electron-free zone after about 
one plasma oscillation, and eventually accelerated by the huge electric field 
produced by the positive ion background with gradients up to TV/m. Up- 
to-date experimental verifications demonstrated gradients of the order of 
100 GeV/m, that can potentially accelerate electrons up to energies in the 
GeV-range within a few millimeters. 

In this paper we study the important issue of longitudinal wake fields 
produced within electron beams accelerated with high-gradient fields. We 
assume that the acceleration distance d a is much smaller than the the over- 
taking length. This is the distance travelled by the electrons as a light signal 
from the tail of the bunch overtakes the head of the bunch. Given a bunch 
of rms length o z , the overtaking length can thus be written as 2y 2 o z , and 
corresponds to the radiation formation length 2y 2 A calculated at A = o z , 
A = A/(2n) being the reduced radiation wavelength. When d a «: 2y 2 o z , the 
electrons can be assumed to be accelerated at a single position z A down 
the beamline. This is the case for laser-plasma devices, since acceleration 
in the GeV range takes place within a few millimeters only. However, it 
is not the case for conventional accelerators, that feature typical gradients 
up to a few tens of MeV/m, and thus need several tens of meters to reach 
the GeV range. The assumption d a «: 2y 2 o z greatly simplifies wake calcu- 
lations. In particular, when this condition is verified, the wake generated 
along the part of the trajectory following the acceleration point z A is inde- 
pendent of any detail of the particular realization of the accelerator. In this 
sense, our study is fundamental, because it remains valid independently of 
the particle accelerator technology chosen, provided that d a «: 2y 2 o z . One 
may also have contributions to the wake generated along the part of the 
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trajectory following the acceleration point Za- These contributions depend 
on the physical nature of the accelerator, can be separately calculated, and 
will be neglected in this paper because they do not affect the bunch in the 
case of a laser-plasma accelerator. 

We base our study on the solution of Maxwell's equations for the longitudi- 
nal field. With the help of the fundamental example of an electron travelling 
in uniform motion we show that the paraxial approximation can be ap- 
plied to describe electromagnetic sources up to the observation point, while 
sources after the observation point in the beam propagation direction can be 
neglected. Then we make use of the paraxial approximation to calculate the 
Fourier transform of the longitudinal field produced by fixed electromag- 
netic sources (current and charge densities) at a certain observation plane 
down the beamline. This expression is at the basis of our treatment, because 
it allows to calculate a closed expression for the longitudinal impedance and 
for the wake function to be evaluated numerically. The knowledge of the real 
part of the impedance further allows verification of the energy conservation 
principle. Impedance and wake function yield an analytical expression in 
the asymptotic limit of a large distance of the electron beam from the accel- 
erator compared with the overtaking length, i.e. in the far-field zone for all 
wavelengths of interest (up to A ~ a z ). This asymptotic limit is of practical 
relevance, as it allows simple estimations of the impact of the longitudinal 
wake on the electron beam energy change. 

We give an application of our theoretical work by studying the feasibility 
of Free-electron Lasers (FELs) drivers, proposed as a possible use of laser- 
plasma accelerators. 

Coherent sources of electromagnetic radiation have become very important 
research tools in science and industry. Recent advance in conventional (RF) 
particle accelerator techniques allow to construct and operate Free-Electron 
Lasers in the VUV as well as in the X-ray range. Lasing at wavelengths 
shorter than the ultra-violet can be achieved with a single-pass, high-gain 
FEL amplifier. Because of the lack of powerful, coherent seeding sources, 
short-wavelength FELs work in the so-called Self- Amplified Spontaneous 
Emission (SASE) mode, where the amplification process starts from shot 
noise in the electron beam Experimental realization of SASE FELs 

developed very rapidly during the last decade. The shortest wavelength 
ever generated by an FEL, A = 13 nm, has been achieved in 2006 at FLASH 
(Free-electron LAser in Hamburg). Regular user operation of FLASH started 
in 2005 [7]. Currently, this facility produces GW-level laser-like radiation 
pulses with 10 to 50 fs duration in the wavelength range from 13 to 45 nm. 
Recently, the German government, encouraged by these results, approved 
funding a hard X-ray SASE FEL (XFEL) user facility, the European XFEL ||8|]. 
The US department of Energy (DOE) has approved the start of construction 
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of the Linac Coherent Light Source (LCLS) at the Stanford Linear Accelerator 
Center (SLAC) [9]. The LCLS and the European XFEL project are scheduled 
to start operation in 2009 and 2013 respectively. 

With the recent progress in the development of laser-plasma accelerators 
lUi 0r 1st] a discussion started about wether and how this technology can be 
used to provide a cost-effective driver for SASE FELs. The very short length 
scale of plasma acceleration may eventually allow for the development of 
very compact laser-like sources both in the VUV and in the X-ray range. In 
this paper we will not address the delicate issue of technical realization of 
such table-top FEL drivers. These should finally provide short bunches up 
to about 10 fs duration in the 100 kA current class (i.e. with a charge of about 
1 nC), featuring a relative energy spread of 0.1% for beam energies in the 
GeV range and a normalized emittance in the order of 1 mm-mrad. Instead, 
we will restrict our discussion to fundamental electrodynamical questions, 
supposing that laser-plasma drivers already exist. Once this assumption is 
made, there are two ways of taking advantage of such device. 

The first way is conservative. Up-to-date laser-plasma technology can pro- 
duce beam energies up to about 1 GeV. With this energy, one may use 
parameters for the system electron beam-undulator in the same range of 
those at FLASH. In this fashion one may eventually achieve wavelengths 
up to 5 4- 10 nm. This application amounts to substitution of the conven- 
tional accelerator system at FLASH with the table-top laser-plasma driver. 
This would not lead to a table-top VUV-FEL, but it would certainly reduce 
its size. The only condition for this scheme to work is that a laser-plasma 
driver capable of providing the same beam parameters as the accelerator 
system at FLASH must first be built in reality. 

The second way is more ambitious, and consists in using laser-plasma accel- 
erators in the GeV range to build Table-Top X-ray FELs by using a completely 
different set of parameters for the system electron beam-undulator. Namely, 
an order of magnitude smaller electron energy, an order of magnitude larger 
peak current, an order of magnitude smaller undulator parameter and an 
order of magnitude shorter undulator period. The present study allows to 
estimate the influence of longitudinal wake fields on the relative energy 
change of the electron beam. In particular, we estimate the magnitude of 
the induced correlated energy-change as a function of the position within 
the bunch, as electrons travel through the undulator. This is a fundamental 
effect that cannot be avoided by fine tuning of the parameters of the setup. 
Moreover, the energy change grows during the passage of the beam through 
the undulator, so that one cannot assume a fixed correlated energy change 
through the undulator. Our conclusion is that the magnitude of this effect 
poses a serious threat to the feasibility of a table-top FEL in the x-ray range 
(TT-XFEL), as well as in the VUV range (TT-VUV FEL). 
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Fig. 1. Geometry. 

Our work is organized as follows. First, in Section 2 we calculate the lon- 
gitudinal field in the space-frequency domain. In the following Section 3 
we introduce the basic quantities impedance and wake function, and we 
review their relation with the Poynting theorem. In Section 4 we present an 
analytical expression for the asymptotes of the longitudinal impedance in 
the limit for a large longitudinal distance of the beam with respect to the 
overtaking length. We further verify the energy conservation principle in 
Section 5. Section 6 is dedicated to the calculation of an analytical expression 
for the asymptotes of the longitudinal wake, always in the limit for a large 
longitudinal distance of the beam with respect to the overtaking length. 
Finally, our study about the feasibility of the TT-VUV FEL and the TT-XFEL 
is given in Section 7. Section 8 concludes our work with a few final remarks. 



2 Calculation of the longitudinal field in the space-frequency domain 



Consider the scheme in Fig. 1, representing the object of our study. Consider 
position z A , representing the exit of the plasma accelerator. At this position, 
the nominal Lorentz factor of electrons is y. The accelerator switches on 
both harmonics of the sources and of the field. Since we consider the wake 
generated along the part of the trajectory following the acceleration point 
Za, the nature of the switcher is not important. It is important, though, that 
the switching process happens on a longitudinal scale d a «: 2y 2 A, repre- 
senting the distance that a light signal has to travel before it overtakes an 
electron of a distance A = c/co, c being the speed of light in vacuum and co 
the angular frequency of the radiation in the space-frequency domain. Sig- 
nificant emission of radiation in vacuum is present for wavelengths longer 
than the longitudinal bunch length a z , up to transverse beam sizes a ± < yo z 
as we will see. As a result, for typical ultra-relativistic beams the condition 
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d a «: 2y 2 A can be read d a <§: 2y 2 o z , 2y 2 o z being the overtaking length. 

We begin our investigation with the help of general knowledge of electro- 
magnetic theory | |lo)l , bycalculating the longitudinal electric field in the 
space-frequency domam_| E z (z ,r Xo ,co) produced at frequency co by given 
electromagnetic sources in vacuum and detected at longitudinal position 
z > Za and transverse position r (see Fig. 1). 

The field in the space-frequency domain E(z , r± , co) is the Fourier transform 
of the field in the space-time domain, E(z ,? Xo ,t). As is well-known, for a 

fixed frequency co = co the field in the space-frequency domain is strictly 

■± 

related to the amplitude E a , D (z ,r Xo ) of a monochromatic field oscillating in 
time with angular frequency co. In fact 

_> -± 

E(z , r Xo , t) = E Wo (Za, j±o) exp [-ico t] + C.C. , (1) 

where "C.C." indicates the complex conjugate of the preceding term. It 
follows that 



E(z , r ±0 , co ) = 2nE w (z , r ±0 )6(co - co ) . (2) 

Therefore, one may use the field in the space-frequency domain for a given 
frequency and yet think of the amplitude of a monochromatic field or vicev- 



1 In this paper we define Fourier transform and inverse transform pair of a function 
f{t) as 

f(co) = J dt f{t) exp [icot] ; fit) = J dco /(cu) exp [-icot] . 

For future use, we also write, explicitly, the definitions of the two-dimensional 
Fourier transform and inverse transform of a function gir) in agreement with the 
one-dimensional notation: 

g$) = J drgir)exp[ir-k\; gir) = ^ J dk g(£)exp [-ir ■ k\ , 

the integration being understood over the entire plane. If g is circular symmetric 
we can introduce the Fourier-Bessel transform and inverse transform pair: 

CO oo 

£(/c) = 2tz J drrgir)J ikr); gir) = ±- J dk kg(k)J (kr) , 



r and k indicating the modulus of the vectors r and k respectively, and J being the 
zero-th order Bessel function of the first kind. 
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ersa. In fact, although they have different dimensions, the knowledge of one 

quantity fully identifies the other through Eq. (2). As an example, E satisfies 
the Helmholtz equation 

c 2 V 2 E + co 2 E = Anc 2 Vp-AnicoJ r (3) 

-± 

where p(r, co) and j(r, co) are the Fourier transforms of the space-time domain 
charge density p(r, t) and current density j(r, t). In the same way, an equation 

for E a , is found substituting ; and p with j a , and p w , related to ; and p through 
an analogous of Eq. (2). 

Eq. (3) can be solved with the help of an appropriate Green's function 
G(z - z', r ±0 - r\) yielding 



oo 

E z (z 0/ r ±0 , co) = -4tt J dz' J dr'^ |^/ 2 - ^ j G(z - z' , r ±0 - f ± ) , 



(4) 



the integration in dr' ± being performed over the entire transverse plane. An 
explicit expression for the Green's function to be used in Eq. (4) is given by 



exp < i(co/c) 



G(z - z', r ±0 - r\) = ■ 



\r± ~ r'J\ + (z -z') 2 



1/2) 



4n 



\r± -r'A +(z -z') 



1/2 



(5) 



that automatically includes the proper boundary conditions at infinity. 

Since electrons are moving along the z-axis, we write the harmonic compo- 
nents of the charge and current density, p and j z as 



p (z! , r\, co) = p (r\) f(co) exp 



.coz 

V J 



u(z' - z A ) 



(6) 



and 



J z (z',r' ± ,co) = pcp. (7) 

Notation u(z' - z A ) in Eq. (6) indicates a Heaviside step function centered 
at position z A , whose presence signifies that there are no sources before the 
plasma accelerator, i.e. before the point z' = Za- This describes a switch-on 
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process. When the switching distance d <sc 2y 2 a z , the nature of the switcher 
is not important for the description of the wake associated with Eq. (6) and 
Eq. (7), i.e. the wake generated along the part of the trajectory following 
the acceleration point z A . However, Eq. (6) and Eq. (7) alone violate the 
continuity equation, and should be completed by extra-contributions that 
depend on the nature of the switcher. For example, if one thinks of an 
acceleration process where a low energy bunch with Lorentz factor y is 
accelerated on a distance d up to a Lorentz factor y, one should add to Eq. 
(6) and Eq. (7) the contribution of the harmonic at frequency co associated 
with the beam with Lorentz factor y . Such harmonic is characterized by 
a different longitudinal velocity j6 c with respect to jSc, and by a different 
wave number k = co/v . However, its amplitude is the same of that after the 
acceleration process: 



p (z' f P ± , co) = p (r'±) f(co) exp 



.coz 

1 



[l-u{z'-z A )\ 



(8) 



and 



joz(z',r'±,co) = p cp . (9) 

As a result, the continuity equation is satisfied. However, Eq. (8) and Eq. (9) 
depend on the particular process described (acceleration from v to v), so 
that their contribution depends on the particular switching process selected. 
In this paper we will not consider this contribution. We will focus, instead, 
on the switch-independent part of the problem. Note that for the case of a 
plasma accelerator, Eq. (8) and Eq. (9) should be replaced by contributions 
due to a positive ion current propagating in the opposite direction with 
respect to the electrons. As a result we can neglect such source, that has no 
effect on the longitudinal wake field acting on the electron beam. 

The assumption of separability of variables z' and r\, in Eq. (6), together 
with the fact that the function / is independent of z' may be satisfied as 
long as the transverse electron beam size a ± remains unvaried along a 
formation length. This happens for an angular divergence cr such that 
o ± » 0ey 2 o z . Moreover, the assumption of separability also requires that 
we can neglect dynamical effects of self-interactions, because there is no 
external force acting on the bunch. In this sense, we are developing the zero- 
th order treatment of a perturbation theory, where the perturbation to the 
particles dynamics is given by the self -interaction within the electron bunch. 
The quantity p has the meaning of transverse electron beam distribution; 
in principle, it may depend on the harmonic co, but in the following we will 
assume it does not. Thus, all information about the longitudinal electron 
density distribution f(t) is included in the function f(co), that is its Fourier 
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transform. A typical Gaussian beam model, that will be useful later on, is 
defined by 



2a*. 



(10) 



o ± being the rms beam transverse dimension and by 



/(f) = — z^exp 
y/2na t 



2a\ 



< — > f(co) = (-e)Nexp 



co 2 o 2 



(11) 



where N is the number of electrons in the beam and (—e) the electron beam 
charge. Moreover, o t is the rms bunch duration, connected with the rms 
bunch length by a z = f}co t , so that in terms of lengths 



f(s) = — ^exp 
V2no z 



2ol 



(12) 



We now account for Eq. (6) and Eq. (7) and calculate the derivative of p with 
respect to z' in Eq. (4), but keep the implicit form for G. We obtain 



E 7 = 



4nicof(co) 



y 2 c 



CO 

J dr' ± p (r\) J dz' exp 



.coz 



G(z - z', r ±0 - r' ± ) 

+47T f(co) exp J df ± p (f' ± ) G(z - z A , r ±0 - r?±) . (13) 

The integration domain in dz' in Eq. (13), that is [za, °°), can be represented as 
[za, z ] U [z , oo). The integral in dz' can be written as the sum of two integrals 
in dz' performed over the separate domains \z A , z ] and [z , oo). At this point 
we perform the following operations. First, we neglect the integral over the 
interval [z , oo), i.e. we neglect the effects of electromagnetic sources located 
in [z , oo). Second, we apply the paraxial approximation. This means that we 
solve, for sources located in [z A , z ], the paraxial equation 



c 2 exp 



.CO 
l—Z 

c 



(n 2ico d \ ~ rj -± 
Vi + — — E = 47ic 2 Vp - 4nicoj , 



(14) 



where we introduced the envelope of the field components 
E = Eexp [-icoz/c] , 



(15) 
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because paraxial approximation implies a slowly varying envelope of the 
field with respect to the wavelength A = 2nc/co. Note that the source term of 
Eq. (14) is now multiplied by a phase factor exp[-icoz/c] with respect to that 
of Helmholtz equation, Eq. (3). As a result, we can obtain an expression for 

E z by formally operating in Eq. (13) in the following way. First, we restrict 
the integration limits in dz' to [z A ,z ]. Second, we substitute the Green's 
function G with exp[-icoz/c] ■ G p , where G p is the Green's function for the 
paraxial equation: 



Gp(z z , Tj_ r'x) — 



4tt(z - z') 



7T ex P 



ico- 



• ±0 I _L 



r 



2c(z - z') 



(16) 



We then obtain the following expression for E z : 



E z (z ,r ±0 ,co) = 
icof(co) 



y 2 c 



za 



(z - z A ) 



exp 



10JZ A 

2fc 



J dr' ± p (r'±) 



exp 



icoz' . \r ±0 -r'j_ Y 
2y 2 c lC ° 2c(z - z') 

r± ~ r'± I 2 



ico- 



2c(z - z A ) 



(17) 



We do not present here a general proof of the validity of Eq. (17). We rather 
verify its correctness in the particular case Za — > -°°,z = 0,p o (r\) = &{f'±) 
and f(co) = (-e). This corresponds to Fourier transform of the field of a 
particle in uniform motion, calculated at position (0, rj_ ). In the chosen 
limit, Eq. (17) becomes 



E z (0,r ±o ,eo = — ^ 



u 

! d 4 



exp 



ICOZ 



2y 2 c 



- ICO 



±0 

2cz' 



(18) 



depending on the modulus of r ±0 only. The integral in Eq. (18) can be per- 
formed analytically using -coz' /(2y 2 c) as integration variable in place of z' 
and taking advantage of the following relation, that is valid for values a > 0: 



/f-p['h + f)]-/*^-*(^. 



(19) 



where K n indicates the n-th order modified Bessel function of the second 
kind. We thus obtain the following expression: 
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E z (0, r ±0 , a>) = - — K 

y 2 c \ cy 



(20) 



in perfect agreement (aside for notational differences) with Eq. (13.80) of [10], 
following directly from the Fourier transform of the time-domain electric 
field. 

To conclude this Section, we give an explicit presentation of Eq. (17) in the 
case z A — > -oo, that corresponds to the steady state solution for the electric 
field, and for z A = 0, that represents the switch-on case without loss of 
generality. 

In the case z A — > — oo we can make still take advantage of Eq. (19). Then, in 
the steady state limit, Eq. (17) becomes 



v i -> \ 2la) 
E z (z , r ±0 , oo) = — — exp 

y l c 



icoz 
2y 2 c 



f(oo) J dr\ p (r\)K | 



co\r ±0 - r\ 



yc 



(21) 



It should be noted that in this case, the dependence of E z on z is restricted 
to a phase factor only (steady state solution). In the case z A = 0, direct 
substitution in Eq. (17) yields 



E z (z ,r ±0 ,co) = - 



cof(co) 



y 2 c 



J dr\p (fx) 



( z 

. r dz' 

J (z -z') 



exp 



iooz' . \r ±0 - r' ± 

+ 100- 



2y 2 c 



2c(z - z') 



2 

y z c 
-\ exp 

00Zn 



100- 



2cz n 



(22) 



It can be verified by inspection that, in the limit for z » 2y 2 A, Eq. (22) 
gives back the steady state solution, Eq. (21). Eq. (21) and Eq. (22) or, more 
in general, Eq. (17) can now be used to estimate the impact of the field 
generated by the electron beam on any particle in the beam. We do so with 
the help of the concepts of wake fields and impedances, that are related to 
the Poynting theorem as reviewed in the following Section. 



3 Poynting theorem in the frequency domain, longitudinal wake fields 
and impedances 



In the time domain, the Poynting theorem |10I1 reads, at any time t: 
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J f-EdV + j S-ndA + j t J (w e - w m ) dV 



(23) 



Notations V and A in Eq. (23) respectively indicate any finite volume and 
its surrounding surface, while h denotes the field of unit vectors normal to 
A, directed outwards. In this way, the second integral in Eq. (23) represents 
the flux of the Poynting vector, while the first and the second integrals are 
volume integrals calculated within V. In Eq. (23) we introduced the Poynting 
vector S and the energy densities w e and w m , understanding that they all are 
functions of time, in the following way: 



S = ^(EXB), „ = ^L 



and w m = — 
8n 



B 



(24) 



Eq. (23) is an exact identity following from Maxwell's equations in the time 
domain. Quantities in Eq. (23) have the dimension of [energy ]/[time]. 

As is well-known [10], a version of the Poynting theorem in the space- 
frequency domain can be derived starting with the Poynting theorem in the 
space-time domain formulated for a monochromatic field like in Eq. (1) and 
averaging over a cycle of oscillation in time. Then, use of Eq. (2) yields the 
following complex relation for space-frequency domain quantities: 



2 J J -EdV+ J S-ndA + ico J (w e - w m ) dV = . 



(25) 



Here we define the complex, space-frequency domain version of the Poynt- 
ing vector S and of the electric w e and magnetic zv m energy densities in 



vacuum as 



S = — [ExB 
2n 



xv P = 



471 



and zv m = — 
4n 



B 



(26) 



Eq. (25) is an exact identity following from Maxwell's equations in the 
frequency domain. Note that now quantities in Eq. (25) have the dimension 
of [energy]/[frequency], i.e. of a spectral energy density integrated through 
a given surface. In our case, we can separately write real and imaginary 
parts of Eq. (25) as 



:j Re[;*-E dV + j 



Re 



S ■ h 



dA = 0. 



(27) 
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Cylindrical surface A (section) 




Electron 
beam 



Volume V 



Fig. 2. Geometry for application of Poynting theorem. 



and 



v 



Im 



• E 



dV + co (w e - w m ) dV + Im 



S ■ h 



dA = 



(28) 



Eq. (27) expresses energy conservation of time-averaged quantities, while 

Eq. (28) "relates to the reactive or stored energy and its alternating flow" 

(cited from [jiol]). In our case we consider the volume V and the surface A 

of cylindrical shape as depicted in Fig. 2. From Eq. (27) follows the equality 

between the total energy losses of the whole bunch and the energy of co- 

-± 

herent radiation in the far zone, a fact that we will verify later on. Since j is 
directed longitudinally we can write such equality (always with reference 
to Fig. 2) as 



2 J Re[//E z ] dV=-^-J 



dA , 



(29) 



A, 



where E ± indicates the radiation field. The volume integral of the quantity 
2/ z *E 2 can be explicitly written with the help of Eq. (22) and Eq. (7) as 



2 J j;E z dV = \f(co)\ 2 X(co) 



(30) 
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where 



2 

X{cS) = ~ J dr' ± J dr* ± p '(r* ± )p (r* ± ) J dz' 

o 

r. r dz" 

X \ l J (z'-z") e 



Mz"-/) . |r ±0 -r' ± | 2 
+ zcu- 







7 2 c f. I 
H exp \ ico- 



r_Lo - r'_ 
2cz' 



2y 2 c 



exp 



2c(z' - z") 



zcuz 



(31) 



Note that that p is a real quantity. The quantity X(co) has the property 
X(cu) = X *{—co), as it can be directly verified. This means that Re[X](cu) is an 
even function of co, while Im[X](cu) is odd. X(co) is strictly related to AE to t, 
the total energy lost by the bunch as it travels within the volume V. In fact, 
AEtot is given by 



OO OO 

=^J dA J d <\M r =^J dA J 



AStot = — I dA I dt 

A -oo 

Comparison with Eq. (29) yields 



dco Ej 



(32) 



A -oo 



OO oo 

AStot = J ^ J Re [/;E Z ] dV = ^ J dcv \f(cv)\ 2 X(co) . (33) 



Since |/| 2 is even, Re[X](cu) is even, and Im[X](cu) is odd, we conclude that 
Afi to t is a positive real quantity, as it must be. 

As we will briefly review here, X(co) is also strictly related to the impedance 
of the system. When the bunch is ultra relativistic, the longitudinal impedance 
of the system, Z (co) is typically given as the Fourier transform of the wake 
function G (As), that is 



Z (co) 



oo 



d(As) 



G (As) exp 



ico- 



As 
i — 

Pel 



(34) 



where the wake function is defined as 
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G (As) = -^Ljdz' E z (As, t)\ t=z , m . (35) 

— CO 

Here E 2 (As, t) indicates the longitudinal component of the time-domain elec- 
tric field generated by a source particle acting on a test particle at longitudi- 
nal distance As from the source. This field is integrated along the test particle 
trajectory, and divided by the electron charge (— e), so that e 2 G (As) is the 
energy (gained, or lost) by the test particle due to the action of the source. In 



agreement with [HO we take the test particle behind the source for positive 
values of As. According to the given definition of wake function, one should 
integrate the longitudinal field over the entire trajectory However, there is 
no principle difficulty in considering only part of the trajectory, let us say, 
up to position z. Mathematically, this means that the upper integration limit 
in Eq. (35), i.e. oo, should be substituted with z. In this way, G = G(As, z). 
Moreover, our trajectory is supposed to start at Za = 0, that allows one to 
substitute the lower integration limit in Eq. (35), i.e. -oo, with 0. We thus 
obtain 



z 



G (As,z) = I dz' E z (As,t)\t= z >/(p c ) ■ (36) 



In both Eq. (34) and Eq. (36) we considered a case when test and source 
particles move along the longitudinal z axis. More in general, we should 
include given transverse offsets of test and source with respect to the z axis. 
Therefore we have to modify Eq. (34) and Eq. (36) to include a dependence 
on such offset, i.e. on the test and source transverse coordinates r ± T and r ± s: 
in this way G = G (As, z, r xSl t±t)- In order to make our definitions indepen- 
dent of r ± T and r ± $, we integrate over the transverse particle distribution in 
dr ±T and dr ±s . This can be interpreted as a substitution of test and source 
particles with disks of charge longitudinally separated by a distance As. We 
thus obtain 



G(As,z) = i J dr' x J dr* x p (r' x )p {r 7 ' \)G (As,z,r' x ,r* x ) . 



(37) 



In Eq. (37) we used the fact that p is independent of the longitudinal posi- 
tion. With the redefinition in Eq. (37) we can further consider the impedance 

Z(co,z): 



oo 

Z{co,z) = J — — G(As,z)exp 



. As 
ico— 

fie 



(38) 
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As is well-known, since G(As, z) is a real function one has Z(co, z) = Z *(-co, z), 
i.e. the impedance has even real part and odd imaginary part. Thus, if 
we split the wake function in the sum G = G s + G A of a symmetric and 
antisymmetric part, the Fourier transform of G$ gives back Re[Z], while the 
Fourier Transform of Ga yields Im[Z]. The symmetric part of the wake, Gs, 
is also called the active part and is related with the energy lost by the bunch 
through radiation. The antisymmetric part Ga instead, is called the reactive 
part and is related with energy redistribution within the bunch, but not with 
the energy radiated. This can be seen by writing the total energy lost by the 
bunch in terms of G(t), where T = s/(/3c): 



oo oo 



A8 tot (z) = Jdtjdz G(T,z)/(f)/(f-x) = j dz G(t,z)^[/](t), 



(39) 



-OO —00 



where J?1[/](t) indicates the autocorrelation function of /. Since ^[/](t) is 
even in t, only the symmetric part Gs enters effectively in the expression for 
Afitot- Finally, using the relation 



j dtgi(t)g 2 *(t) = dco gl ((o)g*((v) (40) 

— oo — oo 

to simplify Eq. (39) we can write 

oo 

A£ tot (z) = ^fdco Z(co,z) \f(cof , (41) 



where we used the autocorrelation theorem to find jfl[f] = \f\ 2 = ( ^[/] ) * 
and Eq. (38). In the case |/(cu)| 2 = e 2 (single electron), by definition of 
impedance, Eq. (34), we have 



Z(co,z) = - J z *E z dV (42) 



v 



and from Eq. (30) we have 



Z(co,z) = ^X(co,z), (43) 
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that yields a practical algorithm to compute both real and imaginary parts 
of the impedance (and, subsequently of the wake), through the evaluation 
of Eq. (31). 

Finally, a word of caution. The impedance formalism is often used in order to 
describe coupling impedances. A typical situation is the following. A given 
electron in a bunch induces charge and current densities in the vacuum 
chamber walls. In turn, these charge and current densities are responsible 
for spurious electromagnetic fields acting back on electrons following the 
exciting one, thus generating a coupling wake function (and thus its re- 
lated coupling impedance). This coupling wake function G c must vanish 
in front of an exciting particle in the limit for y — > oo, i.e. for electron ve- 
locities v — > c. As a result one can take G c (As) = for As < (see 
This property is necessary (although not sufficient) to demonstrate that the 
coupling impedance Z c (i.e. the Fourier transform of G c ), obeys the well- 
known Kramers-Kronig relation 1 12|, 13], thus providing a link between real 



and imaginary parts of the impedance. However, our impedance does not 
obey Kramers-Kronig relations. In fact, we are dealing with a mechanism 
of generation of the wake function that is fundamentally different from that 
of coupling wakes. In particular G(As) can be different from zero both for 
As < and As > 0, a result of the fact that we are dealing both with tail-head 
as well as with head-tail interactions in the time domair{H 



4 Analytical asymptotes of the longitudinal impedance 



Let us write an explicit expression for the impedance with the help of Eq. 
(31): 



Z(co, z) = Z x {co, z) + Z 2 (co, z) . (44) 
Here we have defined, for computational convenience, the quantities 



Z 1 (a>,z)=~ J dr' ± J dT*j_p;(i $ Jp (T* ± ) 

C J , r dz" \ico{z" - z') 

X dz - -exp 

J J {z'-z") 







2y 2 c 



+ ico- 



r' i - f' 



2c(z' - z") 



(45) 



2 Note that causality with respect to time remains obviously valid but it cannot 
be exploited, as in the case of coupling wakes, giving a relation between real and 
imaginary part of the impedance. 
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and 



Z 2 (co,z)=-J dr' ± J df'^ir'Jpoir"^ 

z 

x I dz' — 

J z' 



. \r\-r'\ 





icoz' 


exp 


2y 2 c 



(46) 



Aside for particular cases, the task of calculating real and imaginary part 
of the impedance must be performed numerically. It is possible, however, 
to obtain important analytical information about the impedance and wake 
function in the asymptotic case z » 2y 2 A from Eq. (44). 



4.1 Real part 



Let us consider Re[Z] for z » 2y 2 K. We first deal with Z 2 . With the help of 
Eq. (19), it can be shown that in the limit for z » 2y 2 A, Z 2 tends to a real 
quantity independent of z: 



Z 2 (co) 



= -2cJ d r\J 



dr" ±Po *(r' ± )p (r" ± )K 



co | r'± - r' 



(47) 



The integral in df'x_ can be interpreted as a convolution product between 
two functions of f ± , i.e. p and K . This is equal to the anti-Fourier transform 
in two dimensions of the product of the spatial Fourier transforms of p 

and K , that we will call respectively p (k) and K (k), k being the conjugate 
variable to r ± . By exchanging the integration in dk and that in dr' ± one 
obtains the following alternative representation of Z 2 : 



Z 2 (q>)=~ J dk\p (k)\ 2 K (k) . 



(48) 



As concerns Z x we make use of the two identities 



Ink z'z" 
CO z' -z' 

^ dr ± exp 



exp 



ico\r' ± - r" 



2c(z' - z") 
ico\r ± - fj_\ 2 ico\r ± - r" ± ] 



2cz' 



+ 



2cz r 



(49) 



and 
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Z 2 

J g(z') dz' J g* 



(z") dz" + C.C. = 



z 



(50) 



and of Eq. (19). With the help of these relations we rewrite the real part of 
Zi in the limit for z » 2y 2 /t as an expression independent of z: 



RetZilM = 



or 



J* / 



df±p (r'±)K 



co \r ± - r'j 

yc 



(51) 



Applying Eq. (40) with gi = g 2 to Eq. (51) and remembering that the spatial 
Fourier transform of the convolution product of two function is equal to 
the product of the spatial Fourier transform of the same function, we can 
rewrite Eq. (51) as 



co 



4n 3 y 2 



J dk\p (k)K (k) 



(52) 



Thus, putting together Eq. (48) and Eq. (52) we obtain 



Re[Z]M = J dk\p (k)K (k}\ 2 -^ j dk\p (k)\ 2 K (k) . (53) 

Substitution in Eq. (53) of the explicit expression for K (k) 



K (k) = J dr ± exp\ik-r ± ]K 
2nc 2 y 2 1 



/ -* i \ 
co r i 



yc j 



CO 

I 



= 271 dr ± r ± K 



cor ± 



\Jo(kr ± ) 



[l + (cyk/co) 2 ] ' 



(54) 



together with the new notation 6 - ck/co yields for the sum of two integrals 
in Eq. (53): 



Re[Z](<u) 



-'J 



dd \p o {0,co) 



y 4 d 2 



(1 + y 2 2 ) 



2 " 



(55) 



Up to now we did not make use of any particular model for the electron 
beam. Choosing Eq. (10) as a model for p we can substitute into Eq. (55) the 
following expression for p : 
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Fig. 3. Plot of the universal function C, Re , as a function of a. 



p (O, <y) = - exp 



d 2 M 



2c 2 



This gives 

Re[Z] = Z C Re (a) , 
where 



CRe(«) = {l - (1 + a 2 )exp[a 2 ]r(0,a 2 )} , 



the dimensionless parameter a is defined as 



a = 



CO0j_ 

yc 



and 
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oo 

r(0,x) = j jexp[-f] 



(60) 



and Z = An/c indicates the free-space impedance. A plot of the universal 
function C,R e (a), that is the real part of the impedance in units of Z 0/ is given 
in Fig. 3. Note that Re[Z] exhibits a singular behavior for o±_ — > (i.e. for 
a — > in Fig. 3). This is linked with our particular model, where we chose 



d a -^0. 



4.2 Imaginary part 



Since, in the limit for z » y 2 A, Z 2 is a real quantity, from Eq. (45) follows 
that the imaginary part of Z is 

Im[Z](o;,z) =-^2 Im {/ dr ^/ ^" ±P ;(r\)p (r" ± ) 

f ' , f ' „ 2i \ico(z"-z') . \f\-r 7 \\ 2 

x \ dz' \ dz"- rexp — -—■ + ico-— f- 

J J (z' -z") r I 2y 2 c 2c(z' - z") 







(61) 



In the asymptotic limit for z » 2y 2 A one finds the same result that one would 
have found calculating the impedance with the steady state expression for 
the field, Eq. (21), and substituting z with zero that is 

Im[Z](o;,z) ImjJ df ± J rfr" ±P ;(r\)p (r"J 

Z oo 

f , , f , „ 2i \ico(z"-z') . |' 

x dz dz"- rexp — ^ — + ico- 

J J (z'-z") F [ 2y 2 c 







2c(z' - z") 



(62) 



With the help of Eq. (19) we find the analogous of Eq. (51) for the imaginary 
part of the impedance in the asymptotic limit for z » 2y 2 K: 



lm[Z](cv,z) = -j^ J dP ± J dr'\p;(f 1 _)p (r'' 1 _)K ( ^ 1 ^- ^ ' j (63) 

Similarly as for Eq. (47), we observe that the integral in r" ± can be interpreted 
as a convolution. In analogy with Eq. (48) we obtain 
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Fig. 4. Plot of the universal function Cim/ as a function of a. 
lm[Z](co,z) = J dk\p (k)\ 2 K (k) . 

Finally, Eq. (54) together with the notation 6 = ck/co yields 



(64) 



Im[Z](<u) 



1 f 

= — coz I 



dd \p (d,co) 



(1 + y 2 d 2 ) 



(65) 



Similarly as for the real part, use of Eq. (56) allows to give an explicit 
expression for Im[Z]: 



Im[Z] = (Z z) Cim(«) / 
where 



(66) 



\ 

Cim(a) = -^exp[a 2 ] r(0,a 2 ) , 



(67) 



the dimensionless parameter a is defined as in Eq. (59), T was given in 
Eq. (60), Z = 4n/c indicates, as before, the free-space impedance, and z = 
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z/(2y 2 A). A plot of the universal function C,i m (a), that is the imaginary part 
of the impedance in units of Z z is given in Fig. 4. 

Note that one may use different models for the electron beam transverse 
profile to be substituted in Eq. (65). For example, one may consider a stepped 
profile defined by 



Poft) = 4r C Tb (? ± ) , (68) 



nrtc 

b 



where n is now the beam radius and C rb (r±) is a circle function, equal to 
unity for |r J < n and zero otherwise. The analogous of Eq. (56) is 



By substituting Eq. (69) in Eq. (65) and performing integration we recover 
the already known expression for the longitudinal space-charge impedance 
in case of a stepped profile bunch 



4icz 



con 



con 



corf 



1 - —K, — 



cy 



yc ) 



(70) 



in agreement with [14] where an expression per unit length is given, and 
the overall difference of a minus sign fixes the sign of the energy loss (using 
the convention in Eq. (70) it is negative). 



5 Energy conservation 



It is interesting to verify the energy conservation law for the active part of 
the field, that is Eq. (29). 

With the help of Eq. (30), Eq. (43) and Eq. (55) we can immediately write the 
left hand side of Eq. (29) in its final form as 



J Re[j z %] dV=-^\f(co)\ 2 J d6\p (6,co) 



y 4 2 



(1 + y 2 2 Y 



(71) 



We now need to calculate the right hand side of Eq. (29). This is the spec- 
tral energy density of radiation integrated over a transverse plane in the 
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far zone, i.e. at a large longitudinal distance z with respect to the radiation 
formation length 2y 2 A at all wavelengths of interest, i.e. for A > a z . This 
gives z » 2y 2 o z , 2y 2 o z being the already-defined overtaking length. Radi- 
ation is generated because of the acceleration process, and is linked with 
the longitudinal wakes through Eq. (29). Radiation and longitudinal wakes 
should actually be seen as two faces of the same coin. We have already seen 
that major simplifications arise in the calculation of the impedance (and, 
therefore, of the wake) as one assumes a very short acceleration distance d a 
compared with the overtaking length, d a <sc 2y 2 o z . In particular, in this case, 
the impedance does not depend on the details of the acceleration process, 
that can be considered as taking place at a single point. This fact has its 
reflection on the characteristics of radiation, that are also bound to be inde- 
pendent of how acceleration took place, always provided that d a «: 2y 2 A at 
all wavelengths of interest. The acceleration process works as a switch in 
the space-frequency domain, in the sense that it switches on the harmonic 
contents of the electromagnetic sources at a given wavelength of interest. 
As long as it takes place within a distance much smaller than 2y 2 A, such 
switching process may have very different physical realizations. It may be 
due to an ultra-high field gradient, but also to a bending magnet deflect- 
ing an electron beam from an off-axis trajectory to a trajectory along the z 
axis, or to the crossing of an interface between two media with different 
dielectric properties. All these examples produce well-documented types 
of collimated radiation in the long wavelength range (A > a z ). Radiation 
due to longitudinal acceleration is studied in [10], while if the switcher 
is a bend one obtains edge radiation (see |fl5ll among many others, and 
references therein). Finally, passage through an interface between different 
media produces transition radiation. In particular, coherent transition radi- 
ation in the far zone generated at the interface between plasma and vacuum 
in a laser-plasma accelerator (with d a «: 2y 2 A) has first been studied in [16Q. 
Characteristics of these kinds of radiation coincide as long as the switching 



In reference [17] we showed how radiation produced by a single ultrarela- 
tivistic electron can be suggestively described in terms of a laser-like beam. 
Once the waist of this laser beam is specified, the field distribution at any 
position down the z-axis can be found applying Fourier Optics techniques 
(in free-space this amounts to the application of Fourier propagation equa- 
tion). In particular, a single electron created at Za = by the system depicted 



3 From this viewpoint, we do not agree with [16] about the presence of diffraction 
radiation when the electron beam is created inside the plasma. There is a principle 
difference between a case when an electron travels through a foil with finite trans- 
verse size and another when an electron is created at the interface between plasma 
and vacuum. When d a <sc 2y 2 A radiation characteristics should not depend on the 
transverse characteristics of the plasma. 



process obeys d a «: 2y 2 A 0. 
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in Fig. 2 produces a laser-like beam is with a waist located at Za = 0. As we 
have just discussed, since d a «: 2y 2 A, the field distribution at the waist does 
not depend on the realization of the accelerator setup, and coincides with 
that produced by a magnet edge, or by transition radiation. It is given by 
the field distribution associated with a single electron created at Za = and 
travelling along the z axis |17|1 : 



| ±>(?±)= ?iH^^ Kl (^). (72) 
c z y r ± \ cy J 

In our case though, we are considering an electron beam and not a single 

electron. This means that we are not dealing with a single laser-like beam 

but, rather, with a coherent collection of laser-like beams. At the initial point 

z A = each electron is completely characterized by an offset P x and a 

deflection angle cp' with respect to the z axis. These two vectors identify a 

point in the transverse electron phase space. Each electron, characterized by 
-> -> 

the pair (r'j_,cp' ± ), corresponds to a different laser-like beam, whose waist 

-> _> 
is simply tilted of an angle <p' ± and shifted of an offset r' x . Accounting for 

offsets and deflections, Eq. (72) can be written as 



f 2i(-e)co 
Exs(rJ = -^— exp 



.CO 7> (-> -> \ 

i—cp' ■ [r ± -r' x ) 



r, —r 



J r. 



r 



CO \V I 



r 



cy 



(73) 



We now have to average Eq. (73) over the harmonic component of the 

charge density at frequency co, i.e. p in Eq. (6) calculated at z' = 0, at the 

waist position. In principle, p should depend on both cp' ± and r' x . However, 

in our model, it only depends on r' ' ± , since we posed p = p (r'±)- In our 

-» 

situation of interest we can set cp' = in Eq. (73) and obtain the following 
average over p: 



2ico 
c 2 y 



f(co) J dr ±Po (r_ 



co r i 



cy 



(74) 



The right hand side of Eq. (29) in the space-frequency domain, is given 
by calculating the flux of the complex Poynting vector through the surface 
A (shown in Fig. 2). Such flux is equal to the flux of the complex Poynting 
vector associated with the virtual source through the virtual source, because 
the Fresnel propagator conserves the flux of the Poynting vector. Therefore 
the right hand side of Eq. (29) amounts to 
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cu |r_L - r' ± | 

cy 



(75) 



Note that the integral in dr'± is a convolution product. As before, applying 
Eq. (40) with gi = g 2 to Eq. (75) and remembering that the spatial Fourier 
transform of the convolution product of two function is equal to the product 
of the spatial Fourier transform of the same functions, we can rewrite Eq. 
(75) as 



2n j 



E , d A = - 



A 



2af 

ncy 2 47T 2 



dk\p (k)Tik 



(76) 



where 



having used the notation = c/cok. Substituting Eq. (77) into Eq. (76) we 
obtain the final result 



^ J~|fL| 2 dA = -||/M| 2 J dd\p {6,w) 



y 4 d 2 



A„ 



(i + y 2 e 2 y 



(78) 



that coincides with the left hand side in Eq. (71) thus verifying Eq. (29) as it 
must be. 



6 Analytical asymptote of the wake function 



The anti-Fourier transform of the impedance gives back the wake func- 
tion, that in its turn allows one to calculate the energy change per particle 
averaged over the transverse beam coordinates. 

Analytical results can be found starting with the asymptotic limit of the 
impedance for z » y 2 K. 
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6. 1 Symmetric part of the wake 



-ico 



'As" 



The symmetric part of the wake function G$ can be found calculating 

OO 

G s (As) = ^J da) Re[Z](o;)exp 

— oo 

With the help of Eq. (55) and Eq. (56) one can rewrite Eq. (79) as 

oo oo 

y 4 6> 3 

x2 



(79) 



G s (As) 



7ZC J 



d0- 



(1 + y 2 (9 2 )' 











exp 


c 2 


exp 





(80) 



We first calculate the integral in dco, thus obtaining 



G s (As) 



= __L_ U^* 

V^cTi J (1 + y 2 2 )' 



■ exp 



(As) 2 
40 2 cr 2 



Performing the integral in dO and using the notation A£ 
obtains 



G S {AE) = ^H S (A5) 



(81) 

y(As)/a_ L one 
(82) 



where 



H S (A£) = 



i |A£| + 71 



4Vtt I 



1 - 



(A^) 2 



exp 



(A^) 2 



erfc 



,(83) 



and 



erfc(x) = 1 



- I exp[- 
n J 



■t 2 ]dt . 



(84) 



A plot of the universal function H s , that is the symmetric part of the wake 
in units of y/a ± , as a function of A£ is given in Fig. 5. 

Note that Eq. (82) can be verified in the asymptotic limit for (yAs) 2 » cr 2 , 
i.e. (A£) 2 » 1. If we further assume (As) 2 » o\ we obtain the asymptote 
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for two particles separated by a distance As. It may be useful to remind 
that in agreement with 

El, the case of a test particle in front of the source 
corresponds to negative values of As. In this case, the longitudinal electron 
density distribution and its Fourier transform may be written respectively 
as 



fit) = H) 



(85) 



and 



f{w) = H) \ exp 



. As 
2fic 



+ exp 



-ico 



As 
2j6c 



(86) 



Substitution of Eq. (86) and Eq. (56) into Eq. (78) and use of Eq. (32) yields 



AS 



tot 



CO 

" nc J 



dco < 1 + cos 



cu(As) 



CO 

}/ 



ddexp 



eWo 2 



y 4 d 3 



(l + y 2 d 2 ) 2 ' 

(87) 



The first term in parenthesis {...} is related to single-particle radiation in 
the far zone, while the second term in cos(-) may be interpreted as the 
interference term between the field radiated by the two particles in the far 
zone. Let us consider the interference term alone, that will be indicated with 
Afiinf Performing first the integral in dco we obtain 



AS 



2e 2 



int 



J" dd exp 



(As) 2 
4<9 2 ct 2 



y 4 2 



(1 + y 2 6 2 ) 2 ' 



(88) 



In the limit for (yAs) 2 » cr 2 , Eq. (88) can be written as 



Afiint = 



2e 2 



^ dxexp 



x 2 (As) 2 
4o 2 



2c 2 
As ' 



(89) 



where we have performed a change of integration variable to x = 1/6. The 
energy term AS^ is due to the interaction between the two particles, and 
corresponds to the energy dissipated by the system through the active part 
of the wake. 
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Fig. 5. Plot of the universal function H$ as a function of A£. 

Now, from Eq. (81) we have that the energy lost by one particle is e z Gs(As) — > 
-e 2 /(As) when (yAs) 2 » cr 2 . Since the system is composed by two electrons 
and the wake G s is symmetric, the total energy lost by the system is given 
by twice this value. Thus, in the limit (yAs) 2 » cr 2 , the energy dissipated by 
the system due to the active part of the wake is 



2e 2 

2e 2 G s (As) = -— , 
As 

in agreement with Eq. (89), as it must be. 



(90) 



The energy gained or lost by a single particle at position s within the bunch 
due to the active (symmetric) part of the wake, averaged over transverse 
coordinates is given by 



AS s (s) = (-e) j G s (As)/(s-As)rf(As). (91) 

—00 

Note that this expression can also be presented in terms of impedances as 
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A£s(s) 



oo 

2^1 



dco Re[Z](co)f(co) exp 



-ico 



As 



(92) 



which however appears not to yield any computational advantage over Eq. 
(91). Analogous remark holds for the antisymmetric part of the wake. 

With the help of Eq. (12) and Eq. (82), we can explicitly write A£ s /£ , 
S = ym e c 2 being the nominal energy of an electron (with rest mass m e ) as a 
function of £ = ys/o ± : 



Age 
"67 



(0 = 



oo 



d(A£)H s (£-A5)exp 



2rj 2 



(93) 



where we introduced the new parameter r\ = yo z /o ± , and we called the 
Alfven current Ia = e/(mc 3 ) ^ 17 kA and the beam current J max = eNc/( V~27Zff z )- 

From a methodological point of view it is interesting to conclude our treat- 
ment of the symmetric part of the wake with an additional remark. On the 
one hand we have found that electrons in the bunch lose energy under the 
action of radiative interaction forces within an overtaking length ~ 2y 2 o z 
after the accelerator. On the other hand, an observer in the far zone de- 
tecting the radiation pulse may calculate the retarded position of radiators, 
concluding that electrons radiated inside the accelerator. Methodological 
questions of this kind (electrons radiating when the radiative force does not 
work on them and vice versa) are related with other well-known situations 
of interest. For example, they also arise in Coherent Synchrotron Radiation 
(CSR) from bending magnets and were discussed in [18]. 



6.2 Antisymmetric part of the wake 



The antisymmetric part of the wake function G$ can be found calculating 



G A (As) 



CO 

= 2^J 



dco Im[Z](<i>) exp 



-ico I 



'As; 



(94) 



With the help of Eq. (65) and Eq. (56) one can rewrite Eq. (94) as 
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Using the fact that 
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r eWo 2 ± 






[ c 2 \ 


exp 
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d(As) 



CO 

^ dco exp 



one obtains 



eWa 2 
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i s/n c 2 As 
2e 3 o 3 ± 



■exp 
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(96) 



G A (As) = 



z As r 
2 v^i J (T+ 



d0 



y 2 2 )0 : 



exp 



(As) 2 
4<9 2 a 2 



(97) 



Performing the integral in d0 and using, as before, the notation A£ = 

y(As)/aj_ one has 



ynz 

G A {AE,) = !-L Ha (AQ 



(98) 



where 



H A (A£) = -— (AQ \2 
2 V7i 



'iAii 



n exp 



(A^) 2 



erfc 



|A£| 



(99) 



and we redefined z = z/(2y 2 o z ). A plot of the universal function H^, that 
is the symmetric part of the wake in units of yz/a±, as a function of A£ is 
given in Fig. 6. 

Note that Eq. (98) can be verified in the asymptotic limit (yAs) 2 » cr 2 , i.e. 
(A£) 2 » 1. If we further assume (As) 2 » o\ we obtain the asymptote for two 
particles separated by a distance As as before for the symmetric part of the 
wake. As < indicates again the case when the test particle is in front of 
the source. From the Lienard-Wiechert expression for the electric field of a 
single particle moving on a straight line we can calculate the rate of energy 
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Fig. 6. Plot of the universal function Ha as a function of A£. 



change of the front electron 



dS 



dz y 2 (As) 2 



(100) 



Independently, from Eq. (98) we have that H A (AQ — * -2 sign(A£)/(A£) 2 
when (A£) 2 » 1. As a result, at position z the front particle has gained the 
energy 



A&a = e 2 G A (As) = 



2 



y 2 (As) 2 



(101) 



that corresponds to the rate of energy change in Eq. (100) as it must be. 

The energy gained or lost by a single particle at position s within the bunch 
due to the reactive part of the wake, averaged over transverse coordinates 
is given by 
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Fig. 7. Plot of F in Eq. (105) as a function of s/a z for different values of 77. 

00 

AS A (s) = (-e) j G A (As)f(s - As)d(As) . (102) 

— CO 

With the help of Eq. (12) and Eq. (98), we can explicitly write AS A /S , 
S = ym e c 2 being the nominal energy of an electron (with rest mass m e ) as a 
function of £ = ys/a_ L : 



^A { t.) = hfL n Z f rf(A^)H A (^-A^)exp 



(A^) 2 
2r] 2 



(103) 



Note that Eq. (103) is a function of £ but also depends parametrically on r\, 
and may be presented as 



Ab A ls \ I max z r /s \ 
where we indicated the parametric dependence of t] after the semicolon and 
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A plot of Eq. (105) is given as a function ois/o z in Fig. 7 for different values 
of rj. 

We failed to integrate Eq. (93) and Eq. (103) (or Eq. (104)) analytically. How- 
ever, plots for A&s/& and A& A /& as a function of £ can be computed with 
the help of numerical techniques. In the following Section we will give a 
practical example of application of our work. 



7 Impact on the design of a table-top FEL in the VUV and X-ray range 

The foreseen specifications of the next generation laser-plasma accelerators 
are stunning. The high acceleration gradient will be up to TV/m, producing 
very short bunches about 10 fs long in the 100 kA current class (i.e. with a 
charge of about 1 nC). The beam quality is also expected to rival state-of-the 
art conventional acceleration techniques, featuring a relative energy spread 
of 0.1% for beam energies in the GeV range and a normalized emittance in 
the order of 1 mm-mrad. 

One of the obvious and perhaps most exciting applications envisaged for 
these machines is as drivers for table-top FELs, both in the VUV (TT-VUV 
FEL) and in the X-ray range (TT-XFEL) 1191 ]. Estimation of basic scaling pa- 
rameters for Self- Amplified Spontaneous Emission (SASE) FEL applications 
within the one-dimensional ideal case indicate that a system composed by 
a laser-plasma accelerator and a meter-long undulator may undergo SASE 
process in the sub-nanometer range, thus outclassing all existing and fore- 
seen XFEL projects, both in dimensions and costs. 

Unfortunately, the road map towards the practical realization of this ambi- 
tious goal is not problem-free. Let aside other possible problems we focus 
our attention on a fundamental issue that, in our understanding, consti- 
tutes a serious difficulty. Namely, we want to estimate the impact of the 
longitudinal wake on the electron beam quality using the theory developed 
in the previous Sections. This is possible, because the plasma acceleration 
gradient is such that the electron beam to be injected in the undulator can 
be produced within distances d a «: 2y 2 o z , i.e. well within the applicability 
region of our theory. Note that our estimation is completely independent 
of the physical realization of the accelerator (a laser-plasma device, or any 
other physical principle). Within our model, the accelerator is just located 
at a single point at position z A down the beamline. 



(A^) 2 
2r] 2 



(105) 
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Fig. 8. Geometry of the TT-VUV FEL and of the TT-XFEL setups. 

In the following we will consider two sets of parameters informally under 
discussion within the scientific community. The first refers to a test case 
designed to radiate in the VUV range (A = 25 nm) while the second deals 
with the true table-top SASE XFEL, designed to radiate in the hard X-ray 
range (A = 0.25 nm). Both parameter sets rely on the use of a mini-undulator 
with a 3 mm-period, that would allow low electron beam energy (in the order 
100 MeV and 1 GeV for the VUV and the X-ray case). Very high currents, in 
the order of 100 kA, are supposed to be used . The geometry of the system 
is described in Fig. 8. It is important to note that in both study-cases the 
undulator parameter K - 0.5. This means that the longitudinal velocity of 
particles in the undulator is not significatively altered (y 2 = y 2 /(l + K 2 /2) - 
1.125y 2 ). As a result, wake field calculations in vacuum, performed in the 
preceding Sections, can still be used inside the undulator. 

Description of the FEL setups will be limited to the one-dimensional model. 
In this case, the basic scaling parameter is the one-dimensional p m param- 
eter defined as [20] 



2n 2 j K 2 A 2 



-.1/3 



JJ 



(106) 



where j is the beam current density, A w is the undulator period, K is the 
undulator parameter defined as K = eA w H w /(2mn e c 2 ), where H w is the max- 
imum magnetic field produced by the undulator on the z axis. Finally, the 
coupling factor Ajj, for a planar undulator, is given by Ajj = /o(Q) - /i(Q), 
where Q = K 2 /(2 + K 2 ). The main quantities of interest characterizing the 
FEL process can be written in terms of pn> In fact, the one-dimensional 
power gain length of a mono-energetic beam is 



La = 



A, 



471 V3p 



(107) 
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Fig. 9. Relative energy change A£ s /£ as a function of the position inside the bunch 
s/a s in the case of the VUV table top FEL setup. 

while the relative FEL bandwidth at saturation is close to pm and the power 
at saturation is about pm times the electron beam power. 

Let us first deal with the VUV test case, that aims at producing radiation at a 
wavelength A = 25 nm using an electron energy S = 130 MeV and a 0.7 m- 
long undulator with period A w = 3 mm. The electron beam current is about 
60 kA. The longitudinal beam size is taken to be o z ~ 1 fan. The transverse 
beam dimension is changed, instead, from an initial size a± = 1 [im to a final 
size ct_l ~ 30 [im within a few centimeters. In our estimation we will neglect 
the impact of this change of dimension on the longitudinal wake, and take 
± ~ 30 fim from the very beginning. 

The radiation formation length at wavelength A is estimated to be 2y 2 A, and 
our asymptotic calculations for the impedance are valid for z » 2y 2 A. An 
estimation of the wavelengths of interest is given considering the typical 
length scale that enters in the expressions for the impedance. In our case we 
take A — o z = 1 p.vs\, which yields 2y 2 A - 10 cm, a few times smaller than 
the size of the TT-VUV FEL, because the planned length of the undulator is 
about 0.7 m. This allows us to use asymptotic expressions Eq. (93) and Eq. 
(103). The only parameters needed to be plugged into these equations are 
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-0,4 



Fig. 10. Relative energy change A£>a/&o as a function of the position inside the 
bunch s/o s in the case of the VUV table top FEL setup. This plot refers to the 
position z = 0.7 m. 

thus y = 260 and r\ = 8.7 and I = 60 kA. Moreover, since the asymmetric 
part of the wake is proportional to z = z/ (2y 2 a z ), we need to fix a position 
along the longitudinal axis to calculate Eq. (103). Here we set z - 0.7 m, 
corresponding to the foreseen undulator length. The relative energy change 
ASs,a/£>o calculated with the help of Eq. (93) and Eq. (103) is plot respectively 
in Fig. 9 and Fig. 10 as a function of the longitudinal coordinate inside the 
bunch. While travelling down the undulator, a correlated energy change 
develops along the electron beam. In our case the chirp is non-linear, but in 
order to estimate the magnitude of the effect we can use the linear energy 
chirp parameter 



dy 1 

a = —h—r . (108) 



At ycop 2 : 



ID 



The effect of linear energy chirp starts to play a significant role on the FEL 
gain when a > 1. Intuitively, this means that the relative energy change 
becomes comparable with the FEL parameter on the scale of the coherence 
length. Using an estimation of the slope around s = in Fig. 10 and the 
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other problem parameters one finds a 12, that indicates a very large 
effect. It is important to note that in our case the energy chirp in Fig. 10 
depends on the electron beam profile but also on time, because it develops 
along the undulator. This effect is fundamental, it cannot be avoided, and 
is directly linked with the feasibility of the proposed FEL scheme. In other 
words, radiation generated in one part of the undulator cannot interact in 
resonance in another part of the undulator, and the amplification process is 
destroyed. 

In these conditions, our conclusion is that the proposed TT-VUV FEL setup 
will not work. 

Let us now turn to consider a table-top SASE XFEL scenario aimed at pro- 
ducing radiation at a wavelength A = 0.25 nm using an electron energy 
S = 1.2 GeV and a 3m-long undulator with period A w = 3 mm. The elec- 
tron beam current is taken in the 100 kA range. We take again o z = 1 pm, 
and ct_l = 30 ^m. Similarly as before we estimate the radiation formation 
length taking A — o z . This yields 2y 2 A — 10 m, longer than the TT-XFEL 
undulator length, that is about 3 m. This makes it impossible to use our 
asymptotic expressions, Eq. (93) and Eq. (103), because the entire setup is 
well within the formation length. This situation needs further study, based 
on detailed numerical simulations, that goes beyond the scope of this work. 
However, as before, we can say that longitudinal wake fields will be respon- 
sible for an energy chirp that is both profile dependent and time dependent. 
Moreover, even though the electron beam energy has increased of an order 
of magnitude (thus leading to a decrease of the energy change level), the 
Pid parameter is decreased of an order of magnitude. Since the undulator 
length is only a factor three shorter than the formation length, we conclude 
that wake fields constitute a major effect in this case too. Such effect poses a 
serious threat to the operation of the TT-XFEL setup. 



8 Conclusions 



This work constitutes the first study of the impact of longitudinal wake 
fields on the quality of electron beams produced with high field-gradient 
accelerators. We restricted our attention to the analysis of the wake gen- 
erated along the trajectory following the acceleration, and assuming that 
the acceleration happens on a short longitudinal scale compared with the 
overtaking length, i.e. d a «: 2y 2 o z . Thus, our consideration does not depend 
on the particular realization of the accelerator. However, given present tech- 
nological developments, one of the most relevant applications of our study 
should come from the realm of laser-plasma accelerators. 
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We calculated longitudinal symmetric and anti-symmetric parts of the wake 
function as well as real and imaginary part of the impedance with the 
help of paraxial approximation within a space-frequency domain formula- 
tion of Maxwell's equation. While the general expressions for wakes and 
impedances need numerical techniques to be evaluated, it is possible to 
present analytical expressions for the asymptotic limit when the electron 
bunch has reached a position, down the beamline, that is far from the ac- 
celeration point, with respect to the overtaking distance. It should be noted 
that in these cases the wake and the impedance are proportional to universal 
functions, valid for any set of parameters. 

Our results can be used as analytical benchmarks for computer codes. More- 
over, the the asymptotic result for the antisymmetric part of the wake is of 
fundamental importance, in the sense that its validity is not restricted to the 
particular model studied in this paper (fast acceleration). To the best of our 
knowledge, Eq. (98) constitutes the first analytical solution to the problem of 
space-charge wakes (calculated per unit length). Taking advantage of simi- 
larity techniques we presented such result in terms of a dimensionless object 
dependent on the problem parameters (yr]z/o ± ) multiplied by a universal 
function H A (see Eq. (99)). Eq. (98) can further be used to calculate the rela- 
tive energy change of a particle within a given bunch. In case of a Gaussian 
longitudinal profile the relative energy spread reduces to Eq. (104). Plots 
presented in Fig. 7 will help the reader to perform a fast estimation of the 
influence of the space-charge wake at the stage of planning of experiments. 

As a particular example, we used our results to estimate the impact of the 
longitudinal wake fields on the energy change of the electron beam in a 
table top FEL setup, both in the VUV range (A = 25 nm) and in the x-ray 
range (A = 0.25 nm). In both situations, such energy variation is time and 
shape dependent. This effect is fundamental, in the sense that it cannot be 
avoided. The total energy deviation is found to be an important effect, and 
constitute a major problem for the exploitation of the SASE amplification 
process. 
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